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The scattering among three particles interacting through 1/r 2 forces, with opposite charges and 
widely different masses, is studied in a coplanar geometry. The present work shows that at low 
impact velocities the output of the collision presents typical fingerprints of chaos. The details of the 
^\ ' process are investigated. 
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Starting from the works of Gaspard and Rice jjj], chaotic scattering was found to be ubiquitous. 
■ A large amount of work has been done on potential scattering (i.e. scattering of an elementary 

particle on fixed potential centers) Q-^], but also the scattering among at least three interacting 
bodies has received great attention (see the work of Petit and Henon on the scattering between a 
planet and two satellites ||, or the system helium nucleus plus two electrons Q). A monographic 
issue about this subject has been published in the journal CHAOS ^). 

In all these studies the projectile and a target preserve their internal structure during the collision, 
' but one can analyze also reactive collisions. Kovacs and Wiesenfeld |J studied, in a collinear 

geometry, the scattering between an atom and a diatomic molecule: A + BC <-> ABC <-> AB + C. 
More recently, we studied the chaotic behavior in the reaction between a hydrogen atom (proton 
plus electron) and a projectile proton interacting through Coulomb forces JlO| . We analyzed the full 
three-dimensional problem at very low impact energies. We found that the transition from regular 
to chaotic scattering appears when impact velocity v p is reduced to a value below about 1/10 of the 
classical electron velocity v e . 

In this paper we present another investigation on the same subject: we think that it is worth 
consideration since collisions with rearrangement between electrically charged particles are one cur- 
rent topic in atomic physics, both on the theoretical as well the experimental side juj. The phase 
space of a system of three particles moving in 3D is too large to be easily handled; limiting to 2D 
will allow us to do a more detailed and accurate study. We will address the following topics: a) Do 
features of chaotic scattering appear in two-dimensional collisions? b) If so, do they appear in the 
O ■ form of a sharp order-chaos transition, or as a smooth transition? c) Is it possible to detect some 

traces of irregular behaviour also in the heavy particles motion, instead that only when looking at 
the lightest one? d) Finally, an investigation of the dynamics of the electron during the scattering 
is done. It gives an insight on how the discontinuities on the output function appear. 

The final state of the projectile may be defined through a set of parameters {Af} (e.g. angle 
of deflection, final velocity, ...), which are functions of the input quantities {Ai} (e.g. impact 
parameter, impact velocity, . . . ). When the set of values {Af} depends sensitively by {A^}, i.e. a 
finite variation of {Af} is caused by an infinitesimal variation of {^4i}, then the system is chaotic. 

Usually, one studies the so-called scattering functions, which are the output variables as a function 
of only one input variable, keeping the other input variables fixed. If a scattering function is fractal 
then the system is chaotic but the converse is not necessarily true. In fact, Chen et al. Jl2| showed 
that time- independent Hamiltonian systems with more than two degrees of freedom can have chaotic 
sets with nonzero fractal dimension while at the same time the scattering functions do not show any 
fractal property. The scattering function has fractal behaviour only if the Hausdorff dimension D c 
of the chaotic invariant set satisfies the inequality 

D c > 27V - (2 + q) (1) 
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where N is the number of degrees of freedom and q is the number of conserved quantities of the 
system (in our case N = 6 and q = 4). Therefore a positive answer to the above question (a) will 
give us also a lower bound estimate of the Hausdorff dimension of the chaotic set of the system. 
Our scattering problem has the following alternative outcomes: 

H + H + ->H + + H (2) 
-> e + H+ + H+ (3) 
-f H + H + . (4) 

At the end of the collision the electron may be captured by the projectile nucleus (charge transfer), 
may be ionized (ionization), or may remain bound to the target (excitation of the target). The 
values of the masses and charges are m H + = 1836, m e = 1836, Q e = —Qh+ = — 1 (atomic units 
will be used throughout this paper). 

In our geometry the scattering plane is the plane (x,y). The target proton is initially fixed 
at the origin, r< = (0,0), and at rest, v< = (0,0). The position of the electron is r e = 
(r e cos 9 + u> e t, r e sin 8 + uj e i), where r e is the radius, 6 is the initial phase, u> e — v e /r e is the angular 
velocity. Both v e and r e are equal to 1 in our units. The projectile proton is prepared with velocity 
v p in the positive direction of the y-axis and at the position r p = (b, —a), where b is the impact 
parameter and a > must be chosen so large that the target and the projectile may be considered 
initially as non-interacting. In these simulations a ~ 10 was found to be a good choice. 

The total energy of the system is given by 

n 1 2 1 2 1 2 1 1 1 

E = 7; m H+ v t + o m H+ u p + 7; m ^ v c + 1 r ~ 1 r - 1 r • ( 5 ) 

2 2^2 |r f -r p | |r t — r e | \r p - r e \ 

The classical equations of motion are numerically integrated by using a variable-order variable-step 
Adams algorithm (routine NAG D02CJF). 

Please, note that the choice of the outcome is done on the basis of the energy of the electron with 
respect to the two bodies. An algorithm to identify all the possible outcomes is given in ref. 
In particular, in our system three parameters remain free: 8, v p , and 6. Numerical simulations have 
been carried out varying all of them. The output variables {^4/} defined above reduce to just one 
element -the final state of the electron. It is a discrete function of its energy. Discrete variables are 
usually not employed in these studies but are not a novelty: already Bleher et al Jl^] investigated a 
chaotic system with a finite number of possible outcomes. 

In Figures 1 and 2 some consecutive blow-ups of the scattering function are presented as a function 
of b, for two different choices of v p and 6. There are regions where the outcome is a regular function 
of b, alternated with other regions where an irregular behaviour appears. One recognizes a kind 
of self similarity, the same pattern repeating at smaller and smaller scales. This is the distinctive 
aspect of fractal behaviour. The same results appear varying 8, at v p and b fixed. 

A clear picture is obtained by studying the scattering function versus the proton velocity v p . In 
Figure 3 this is shown for one couple of values (b,8). It is visible the great degree of irregularity at 
small v p 's. For higher v p 's the zones where the electron is captured and those where it remains close 
to the target proton are more clear-cut. Ionization is by far the least probable process. One may 
judge that the transition order-chaos does appear quite smooth and the scattering at v ~ 0.2, 0.3 is 
not yet entirely regular. 

In order to quantify the amount of chaos one can compute the fractal dimension of the scattering 
function. Given a reference trajectory with initial impact parameter b, we slightly modify the impact 
parameter to b + e and make another run. The results of the two runs are compared. The starting 
trajectory is said uncertain if the two final states are different. The fraction /(e) of uncertain 
trajectories as a function of the parameter e scales as 

/(e) ~ e 1 -" , for e - , (6) 

where d is the capacity dimension ^| . In Figure 4 we plot the capacity dimension d as a function of 
v p . There is a regular decrease of d for v p < 0.3, then d is nearly constant but positive. 

One of the main purposes of this work is to analyze in detail the behaviour of the electron in 
correspondence of a singularity of the scattering function. In Figure 5 the trajectories of the particles 
are shown for two almost identical runs, differing only by a very small value of the impact parameter. 
We see that initially the electron trajectories in the two cases nearly overlap. The two paths begin 
to differ when the nuclei reach the point of closest approach. The potential felt by the electron may 
be sketched as two wells located around the nuclei divided by a symmetrical ridge. Obviously, the 
location of the ridge dynamically evolves in time as the nuclei move. The electron initially lies in 
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one of the two potential wells and it is able to escape to the other well only if at some time during 
the collision it passes near the ridge with enough velocity to cross it. Note that we are referring to 
the component of the velocity orthogonal to the potential ridge. Once the electron has fallen into 
the other potential well it is very unlikely that it can cross the ridge the opposite way and come 
back to the former nucleus. Two electron differing by an infinitesimal value of the velocity, u, u + e 
will undergo entirely different fates provided that u + e is not enough to cross the potential barrier 
while u is. 

It is interesting to notice that chaotic features appear only as singularities of the scattering 
function: An analysis of the delay time (the time needed by the electron to quit the interaction 
region) reveals that it is not a sensitive function of the scattering. This is related to the fact that 
the electron does closely follow one of the nuclei, and the trajectories of the heavy particles are 
always regular: no complicated patterns appear from a study of their behaviour. 

In summary, a study of the two-dimensional low energy scattering between charged particles has 
been carried out. Simulations have been done using a larger number of runs than in our previous 
paper |^ so diminishing statistical errors. The results confirm the findings of ref. jHJ that the 
scattering becomes irregular when the energy is low enough. An insight on why chaos appears has 
been given by looking at the detailed trajectories of the electron. Abrupt discontinuities in the 
output function are correlated with the structure of the potential seen by the electron. 

It is useful to spend some words about the validity of the model used. At very small impact 
velocities, it is not entirely justified to treat the heavy particles as classical objects, and even less 
the light particle. The processes (2-4) should be studied within the framework of quantum mechanics. 
However: i) one may imagine to apply the same apparatus not to the ground state but to a Rydberg 
state, where using the classical mechanics is justified; ii) the use of classical or semiclassical models 
has recently been extended with satisfactory results in realms thought previously to be amenable 
only to quantum treatment, see for example the classical description of the helium ato m 1151 , or the 
treatment of the hydrogen ionization by electron scattering at energies near threshold |16|T In that 
work, furthermore, the evidence for chaos was found, reminiscent of the results shown in this paper. 
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FIG. 1. Scattering function versus impact^parameter b. 'C means electron capture; T, ionization; 'E', electronic excitation; 
'F', undetermined. Here v = 0.10, 6 = 5.82512 radians. From top to bottom successive enlargements are shown with respect 
to b. The number of runs is about 500 for each plot. 
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FIG. 2. The same as Fig. 1, but now v = 0.120, 9 = 0.25262 radians. The number of runs is about 500 for each plot. 
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FIG. 3. Scattering function versus impact velocity v. b = 0.2893, 
each plot. 
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FIG. 4. Capacity dimension ci versus impact velocity v. Statistical errors are shown. 
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FIG. 5. Upper, trajectories of the particles in the interaction region, v = 0.10, 6 = 5.82512 radians, and b = 2.345. The 
dashed line is the trajectory of the heavy projectile. The electron trajectory is the dotted line. The initial position of the 
electron is close to (0,0). The trajectory of the target nucleus is not shown since it is always close to the origin (0,0). The 
lower panel shows the same process, but now with b = 2.350. 
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